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CO ABSTRACT 
l> 

We apply a Wiener filter to the two-year COBE DMR data. The resulting sky map has 
significantly reduced noise levels compared to the raw data: the most prominent hot and cold 
Q\ spots are significant at the 4-<7 level. Furthermore, the entire covariance matrix of the errors 

in the filtered sky map is known, and it is therefore possible to make constrained realizations 
of the microwave sky with the correct a posteriori probability distribution. The filtered DMR 
sky map is used to make predictions for the Tenerife experiment. Two prominent features are 
O , predicted in a region of the sky not yet analyzed by the Tenerife group. The presence of these 
q features is a robust prediction of the standard cosmological paradigm; if these features are not 

S— < observed, some of our fundamental assumptions must be incorrect. 

Subject headings: cosmic microwave background — methods: statistical 



1 Introduction 



The discovery of anisotropy in the cosmic microwave background (CMB) by the COBE satellite 
(Smoot et al. 1992) marked a major advance in cosmology. CMB anisotropy in general and the 
COBE data in particular now provide some of the most powerful tests of cosmological theories, 
particularly theories concerning the formation of large-scale structure. 

The COBE DMR data have been extensively analyzed, with most of the emphasis being 
placed on the task of estimating the power spectrum of the anisotropy (e.g., Bennett et al. 1994; 
Bond 1995; Bunn, Scott, & White 1995; Bunn & Sugiyama 1995; Gorski et al. 1994; White 
& Bunn 1995) although a few have focused on other problems, such as testing the hypothesis 
that the anisotropy obeys Gaussian statistics (e.g., Hinshaw et al. 1994, Kogut et al. 1995). 
In addition, statistical comparisons have been made between the DMR sky maps and data 
from other experiments (Ganga et al. 1993, Watson & Gutierrez de la Cruz 1993, Lineweaver 
et al. 1995). Less attention has been paid to studying the spatial properties of the sky maps 
themselves to determine, for example, which hot and cold spots on the maps are likely to be 
real, and which are dominated by noise. Such an analysis has the potential to be quite useful 
in comparing the DMR results to other experiments that probe comparable angular scales, such 
as the FIRS (Meyer et al. 1991) and Tenerife (Davies et al. 1992, Hancock et al. 1994, Watson 
et al. 1992) experiments. In this paper we will use Wiener filtering to study these questions, 
extending a similar study that we did for the first-year DMR maps (Bunn et al. 1994). In 
particular, we will use the Wiener-filtered two-year COBE DMR sky maps to make predictions 
for the Tenerife experiment, which probes smaller (although comparable) angular scales than 
COBE. 

There are two major difficulties in studying the DMR maps. The signal-to-noise ratio in each 
pixel is quite low, and pixels near the Galactic plane are contaminated by Galactic emission. 
These problems place strict limits on the accuracy with which the true structure can be recovered 
(Bunn, Hoffman, & Silk 1994; Stark 1993). The techniques of Wiener filtering and "constrained 
realizations" provide one way of mitigating these problems. We assume the correctness of 
the canonical theory of large-scale CMB anisotropy: that the anisotropy forms a Gaussian 
random field with a known power spectrum. We then remove the pixels that are presumed to 
be contaminated, and apply an optimal linear filter to the data in an attempt to clean up the 
noise and see the underlying structure. The technique of constrained realizations helps quantify 
the uncertainties associated with this method. 

The Wiener filter described in this paper is based on techniques developed by Wiener (1949) 
and introduced to astrophysics by Rybicki & Press (1992), and it is similar to a technique applied 
recently to galaxy catalogues (Lahav et al. 1994, Zaroubi et al. 1994) The formalism for making 
constrained realizations of Gaussian random fields was developed by Hoffman & Ribak (1991) 
and has been applied to the study of large-scale structure by Ganon & Hoffman (1993). 

The idea behind the Wiener filter is quite simple. If we are willing to assume that we know 
the statistical properties of the signal and noise contributions to the data, we can apply a linear 
filter that preferentially removes noise and leaves signal. The Wiener filter is the optimal linear 
filter for this purpose. It can be derived in two independent ways. First, it is the optimal 
linear filter in the sense of least squares. That is, the Wiener filter minimizes the expected 
mean-square deviation between the filtered data and the underlying cosmic signal. It is also a 



2 



maximum-likelihood estimator of the cosmic signal. That is, the filtered sky map is the most 
likely realization of the underlying signal at each point. As long as both the signal and the noise 
are Gaussian, the filters chosen by these two criteria are identical. 

Furthermore, we can determine the complete a posteriori probability distribution of the true 
signal given the data and the assumed power spectrum. This probability distribution is Gaussian, 
with mean given by the filtered data. We can compute the covariance matrix of this probability 
distribution and use it to construct "constrained realizations" of the CMB anisotropy with 
the correct a posteriori probability distribution. These realizations are useful for qualitatively 
assessing the uncertainties in the filtered maps, and for using Monte Carlo simulations to make 
predictions for other experiments. 

Although detailed derivations of the Wiener filter have been given elsewhere (Bunn et 
al. 1994, Bunn 1995), we will give a brief sketch of the central results in the next section. We 
will then apply this formalism to the two-year DMR sky maps to identify features in the data 
that have high statistical significance. Finally, we will use the filtered data to make predictions 
for the Tenerife experiment. 

2 Deriving the Wiener Filter 

As usual, we will expand the temperature anisotropy AT in spherical harmonics: 



The signal seen by the COBE DMR experiment is the convolution of AT with a beam pattern 
W. The spherical harmonic coefficients of the DMR signal are therefore W;a; m , where W\ is a 
coefficient of a Legendre polynomial expansion of the beam pattern. The DMR beam pattern 
is often approximated as a Gaussian with a 7° beam width, but we use the more accurate 
characterization given by Wright et al. (1994). 

After removing all pixels within 20° of the Galactic plane, the COBE DMR data contain 
J\f = 4038 pixels. Each pixel contains contributions from both signal and noise: 



Here di is the measured value of the ith pixel, a\ m is a coefficient of a spherical harmonic 
expansion of the actual CMB anisotropy, r,- is the location of the ith pixel on the sky, and n,- 
is the noise. If we denote pairs of indices (l,m) by a single Greek index fj, and define matrices 
Y and W by = Y^(Yi) and W^i = WiS^i [where fj, = (7, to)], then we can write the above 
equation in a more compact form: 



We make the usual assumption of Gaussian statistics: both and n,- are Gaussian random 
variables with zero mean and covariances 







d = Y ■ W • a + n. 




(d-d T ) = c, 

(n ■ n T ) = N, 
(a M n.) = 0. 
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The matrices C and N are diagonal. = C\ is simply the angular power spectrum, which 

can be predicted from a particular theoretical model, and Na = af is the noise variance in the 
ith pixel. [The noise covariance matrix is not strictly diagonal, but the approximation that it is 
diagonal is quite good (Lineweaver et al. 1994).] 

In principle, a and C should be infinite-dimensional, but in practice the COBE beam smooths 
away modes with high /, and we can truncate the spherical harmonic expansion at some / max . 
We shall see below that the filter recovers virtually no information for / ^ 20, and so we generally 
set 'max = 30. The dimension of a is therefore M = (/ max + l) 2 = 961. The covariance matrix 
of the data vector d contains contributions from both signal and noise: 

M = (<f • dF) = Y ■ C ■ Y T + N, (7) 

where the diagonal matrix C = W-C-W contains the beam-smoothed angular power spectrum. 

The Wiener filter is the linear filter F that estimates the true cosmic signal s = Y • a as 
closely as possible, in the sense of least squares. Specifically, if we set 

y = F-d, (8) 

then we want to choose the J\f X Af matrix F to minimize 

A=((y-Y.af). (9) 

Setting dA/dFij = and solving for F, we find the least-squares filter 

F = Y • C • W • Y T -M -1 . (10) 

This derivation makes no use of the fact that the random variables a and n are Gaussian; it 
uses only their covariances. If we make use of our assumption that these quantities are Gaussian, 
we can give an alternative characterization of the filter (10): it gives the maximum-likelihood 
estimator of the signal s. 

Bayes's theorem states that the conditional probability density of a being the correct set of 
spherical harmonic coefficients given the data d is simply f(a\ d) cx f(d \ d)f(a), where f(d \ a) 
is the probability density that the actual data d would have occurred assuming a is the correct 
set of coefficients, and f(a) is the a priori probability of getting the coefficients a. All of these 
probability densities are Gaussian: 

f(a\d) oc exp (~\{d- Y • W • a) T ■ N" 1 • (d- Y • W • a)) exp (-\a T ■ C" 1 • a) . (11) 

By completing the square and performing some unenlightening algebraic manipulations, one can 
show that f(a | d) cx exp ( — \x 2 ), where 

X 2 = (a - < L ) T • (C" 1 + W • Y T • N" 1 • Y • W) • (a - a ML ) + dF ■ M" 1 • d, (12) 

and 

a ML = (CT 1 + W • Y T • N" 1 • Y • W)" 1 • W • Y T • N" 1 • d. (13) 
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It is easy to see from equation (12) that the maximum-likelihood value of a is a ML . The 
maximum-likelihood estimate of the actual CMB signal is therefore s* ML = Y • a ML = F' • d, 
where the filtering matrix is 

F' = Y (CT 1 + W • Y T • N" 1 • Y • W)" 1 • W • Y T • N _1 . (14) 

The two filters F and F' are in fact identical, despite their different appearances. Four our 
purposes, since the number of data points J\f is larger than the number of degrees of freedom 
A4, equation (14) is more convenient computationally, requiring the inversion of a smaller matrix. 

Equation (12) gives us the entire a posteriori probability distribution of a, i.e., the conditional 
probability density of a given the data and our assumed power spectrum. This probability 
density is a multivariate Gaussian with mean a ML and covariance matrix 

K = (CT 1 + W • Y T • N" 1 • Y • W)" 1 . (15) 

This probability distribution tells us how likely a particular value of a is, given our assumptions. 
The square roots of the diagonal elements of this matrix can be interpreted as standard error 
estimates for the elements of a ML , and those of Y • K • Y T as error estimates for the elements of 
s ML = Y-a ML , which is the maximum-likelihood value of the true CMB signal. In addition, armed 
with the full covariance matrix K, we can construct constrained realizations of the microwave sky 
by simply choosing Gaussian random variables from the appropriate distribution. In this way 
we can construct simulated sky maps with a probability distribution that correctly incorporates 
both our presumed model and the existing data. These simulated maps can be used to make 
predictions for other experiments, as we shall see below. 

Unfortunately, equation (3) is not a completely adequate model of the DMR data, because 
the data do not contain information about all of the multipole moments. In particular, the 
data are completely insensitive to the monopole contribution, and the dipole is contaminated by 
the much larger kinematic dipole due to our motion with respect to the CMB rest frame. The 
Wiener filter derived above must be modified to account for this. For simplicity of notation, we 
will ignore beam-smoothing by setting W = 1 throughout the remainder of this section. 

Suppose that the data contain no reliable information about modes with / < /*. To account 
for the monopole and dipole, we should set /* = 1; if we decided to discard the quadrupole as 
well, we would set /* = 2. Then we should modify equation (3) as follows: 

d = Y ■ a + Z • b+ ft, (16) 

where we have excised the elements of a and the columns of Y corresponding to the modes about 
which we have no information. The dimension of a is now M. = (/ max + l) 2 — (/* + l) 2 instead 
(I mas | l) 2 , and the corresponding dimensions of Y and C are reduced in the same way. Z is an 
Af X (/* + l) 2 matrix containing the columns removed from Y, and b is an (/* + l) 2 -dimensional 
vector containing the unknown coefficients of the excluded modes. 

We can still try to choose our filter F to be a least-squares estimator in the sense of equa- 
tion (9). However, we must impose an additional constraint on F: since we have no knowledge at 
all of the true value of 6, we should demand that the result F-d of filtering should be independent 
of b. We can do this by imposing the constraint 

F-Z = 0. (17) 
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If we solve this constrained minimization problem by the usual method of Lagrange multipliers, 
we find that 



F = Y C Y T M _1 (l - Z • (Z T • M" 1 • Z)" 1 • Z T • M" 1 ) = Y • C • Y T • M" 1 • P. (18) 

This equation differs from the original filter (10) only in the projection operator P. This operator 
simply projects the data vector onto a subspace orthogonal to the unwanted modes, where 
orthogonality is determined with respect to the inner product (x, y) = x T ■ M _1 • y. 

The Wiener filter (18) is also a maximum-likelihood filter: the conditional probability dis- 
tribution for a and b given the data is P(a, b \ d) oc exp ( — \x 2 ), with 



The maximum-likelihood filter F' in equation 14 must therefore be modified by including the 
same projection operator P, and the least-squares and maximum-likelihood filters still coincide. 
Furthermore, equation (15) for the covariance matrix K in the a posteriori probability density 
of a is unchanged. 

We conclude this section with a pragmatic note. It is in fact not generally necessary to 
compute the projection operator P in order to properly project out the unwanted modes. Rather, 
one can achieve identical results by simply including the unwanted modes together with the 
others, so that a is once again an (/ max + l) 2 -dimensional vector. We extend the signal covariance 
matrix C by setting = ooS^ in the rows and columns corresponding to the unwanted modes. 
Then the filter F and the covariance matrix K computed in the previous section give precisely 
the right results. 1 

This result is easy to justify heuristically: since we have no information about the coefficients 
6, it is natural to assign them an infinitely wide prior distribution. It is also straightforward to 
verify by a direct computation that the results of the previous section approach the results of 
this section as the unknown elements of C tend to infinity (Rybicki & Press 1992). 

3 Results 

In this section we will apply the filter derived above to the two-year DMR data. Figure 1 shows 
the result of applying the Wiener filter to the two-year DMR data. The data set used to make 
this map is a weighted average of the 53 and 90 GHz ecliptic-projected two-year sky maps. The 
weights were chosen to be inversely proportional to the noise variances in the maps in order 

1 One does not of course set anything equal to infinity in one's code: matrices with infinite entries are notoriously 
ill-conditioned. In practice, if we use equation (f4) to compute F, we need only C -1 , never C, and we can simply 
set the appropriate entries of C -1 to zero. 



X 2 = (a - d ML ) T ■ (C" 1 + Y T • N" 1 • Y) • (a - a ML ) 

+(b- b ML ) T ■ Z T ■ M" 1 • Z • (6 - 6 ML ) + <F ■ M" 1 • I 



(19) 



where 



a ML = C • Y T • M" 1 • P • d, 
&ml = (1-P)-<£ 



(20) 
(21) 
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to minimize the noise in the average. The usual 20° Galactic cut was imposed on the map 
before filtering. The input power spectrum was a standard cold dark matter power spectrum 
with normalization Q = 20 /iK. Varying the input power spectrum leads to slight variations in 
the amount of small-scale power in the filtered map, but the large-scale features are essentially 
unchanged. It comes as no surprise that there is less information in the Galactic plane than at 
high Galactic latitude. 

The root-mean-square anisotropy level in the filtered map in Figure 1 is 30.6 /iK. The 
ensemble- average expected signal for the power spectrum we have chosen is 42.5 /iK, so we have 
reason to believe that we have recovered a significant fraction of the power. We can assess the 
uncertainty in this filtered map in a variety of ways. Given the data and the assumed power 
spectrum, the a posteriori probability distribution of the true signal is Gaussian with mean 
equal to the filtered signal. The simplest way to assign uncertainties is simply to compute the 
standard deviation of this probability distribution. In other words, we compute the square roots 
of the diagonal elements of the covariance matrix Y • K • Y T . The ratio of the signal to the 
uncertainty is shown in Figure 2. 

We can get a qualitative idea of the statistical significance of features in the map from 
this figure: individual pixels in the most prominent hot and cold spots are significant at the 
three to four sigma level. Of course, the various pixels are not independent: the covariance 
matrix Y • K • Y T has off-diagonal elements. The correlation between the errors in nearby 
pixels is approximately Gaussian with standard deviation 5?7. We can also get a handle on the 
degree of uncertainty in the filtered maps by making constrained realizations of the temperature 
anisotropy in the manner described in the previous section. Four such realizations are shown in 
Figure 3. 

We can define the angular power spectrum of the filtered map: 



This power spectrum, shown in Figure 4, s a useful indicator of the amount of information 
recovered by the filtering process. On large scales the recovered power A\ is comparable to the 
input power spectrum, while for modes with high / the power is nearly zero. The Wiener filter, 
by construction, returns zero values for modes for which we have no information: after all, in 
the absence of any information from the data, the a posteriori probability distribution for a\ m 
is the same as the prior distribution, which is a Gaussian of zero mean. 

4 Predictions Based on the Wiener-Filtered Map 

One use for the Wiener-filtered DMR maps is to make predictions for other experiments. We 
can process the filtered map through the window function of the other experiment to determine 
the most probable signal for the experiment to observe, given the DMR data and the assumed 
power spectrum. Furthermore, since we know the entire conditional probability distribution for 
AT given the DMR data and the assumed prior distribution, we can compute the probability 
distribution for the predicted signal. Since we have incorporated the DMR information, this 
distribution will in general be narrower than the a priori probability distribution associated with 
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our theoretical model. We should therefore be able to place tighter constraints on our model than 
we could without using the DMR data. By using the filtered data rather than the raw data, we 
reduce the uncertainty in the prediction as well as avoiding the technical difficulties associated 
with comparing data sets of different resolutions and beam-switching strategies (Lineweaver et 
al. 1995). If one has access to both data sets, of course, then both can be used as inputs to the 
Wiener filter, rather than using the data from one experiment to make predictions for the other. 

Of course, we only expect to gain significantly if we choose an experiment whose window 
function has significant overlap with the DMR window function, so that the modes that we have 
estimated well from the DMR data make a substantial contribution to the signal in the other 
experiment. The Tenerife experiment would seem to be ideal for this purpose, since its beam 
size of 5°2 (FWHM) is only slightly smaller than the 7° DMR beam. Unfortunately, the fact 
that Tenerife is a triple-beam experiment works against us, since the sensitivity to modes with 
low \m\ is quadratically suppressed. This fact significantly reduces the amount of information 
we can recover from the filtered DMR map. 

In Figure 5 we show the COBE and Tenerife window functions Wi- Wi is the window function 
used to compute the mean-square signal expected in an experiment. Specifically, if i?(r) is the 
response of the experiment when pointed at a particular position r, then 

(R 2 ) = ^J2(2l+l)CiW l . (23) 

Note that this is different from the quantity W\ that we defined above: W\ is used to compute 
the response i?, not the mean-square response (i? 2 ). For COBE, the two are related in a simple 
manner: Wi = Wf . For Tenerife, the beam-chopping breaks the circular symmetry, and so W 
depends on both / and m (see below). In this case, 

Air ' 

w '= 2TTI ^ w ^ Ylm ^ 2 - < 24 ) 

m = — l 

The Tenerife experiment has taken data covering the range 140° to 260° in right ascension 
and 35° to 45° in declination. Qualitative maps of the entire data set have been published 
(Rebolo et al. 1995), but only the strip near 40° has been quantitatively analyzed (Watson et 
al. 1992). We will attempt to predict as much as we can of the Tenerife signal over the entire 
range. The strip of sky observed by the experiment is shown in Figure 1. 

The process of making predictions is fairly simple. We take the spherical harmonic coefficient 
a ML and multiply it by the Tenerife window function W^ Te11 ^: 

WJ en) = (1 " cosma)exp(-i4 en /(/ + 1)). (25) 

Here oren = 2°2 is the Tenerife beam size, and a sin 9 = 8°1 is the amplitude of the beam chop. 2 
The resulting vector contains the coefficients of a spherical harmonic expansion of the most 

2 The window function is only simple in a coordinate system in which the chop is in the azimuthal direction. 
We must either compute a ML in these coordinates in the first place or else rotate to them afterwards. The Tenerife 
experiment always chops in right ascension, so the rotation is easy to perform. 
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probable a posteriori Tenerife signal Rml- 



(26) 



The result is shown in Figure 6. The portion of sky plotted corresponds to the region of 
sky covered by the Tenerife experiment. The root-mean-square value of Rml over this region is 
16.3 //K. For the standard CDM model we are considering, the expected r.m.s. signal is 32.7 /iK, 
and so it appears that we have predicted a significant fraction of the total signal. Recall that the 
Wiener filter returns low values for modes for which it has little information, and values near 
the expected value for modes with high signal-to-noise ratio. The ratio of the Wiener-filtered 
prediction to the a priori expectation value is thus a rough indicator of how much signal has 
been recovered. 

We can of course do much better than this qualitative argument. We know that the vector a 
is Gaussian with mean a ML and covariances given by (15). The Tenerife response R is therefore 
Gaussian distributed with mean Rml an d covariance 



The standard error associated with each point in Figure 6 is therefore <r(r) = ^/Cov(i?(r), R(r)). 
The standard error ranges from 26.4 //K to 29.7 //K. 

One way to get an idea of the amount of information we have gained by this exercise is to 
compare these standard errors with the level of uncertainty we would have without the DMR 
data. If we made no use of the DMR data, then our probability distribution for R{r) would 
simply be Gaussian with zero mean and standard error given by the r.m.s. expected signal 
of 32.7 /iK. Using the DMR data reduces the standard error by about 15%. This amounts to 
at least a modest improvement in our ability to assess the consistency of this model with the 
Tenerife data. Furthermore, some particular features in the prediction, particularly the two 
positive features at large right ascension, have quite high significance. 

In Figure 7 we have plotted the predicted signal and associated errors in several individual 
strips for greater ease of comparison with the actual data. The maximum at right ascension of 
250° coincides almost exactly with the largest peak in the Tenerife map (Rebolo et al. 1995). 

The standard errors simply represent the diagonal part of the covariances in equation (27). 
For a complete statistical description of the prediction, we need to know the correlations between 
distinct points as well. Because of the beam-switching strategy, these correlations can be quite 
important. The general shape of the correlations is as one would expect from the beam-switching 
patterns: a positive lobe flanked by two symmetric negative lobes. The following fitting formula 
reproduces the three main lobes to about 5% accuracy, although it misses small ripples at high 
separation. 



where AS is the declination difference between ri and r*2, and A A = A^cos^i is the separation 
in azimuth. 



Cov(R(h),R(h)) = J2 w F en %(^) w l Ten)Y »(^) K i 



(27) 




(28) 
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As with the original Wiener-filtered map, it may be helpful to look at constrained realizations 
of the statistical fluctuations about the maximum-likelihood Tenerife prediction in order to get 
some idea of the level of uncertainty. Five such realizations are plotted in Figure 8. We show 
both greyscale maps for comparison with Figure 6 and linear plots of the strip at declination 
35° for comparison with Figure 7. 

The Wiener filter depends on the input power spectrum assumed. It is therefore of interest 
to ask how robust the features in the Tenerife predictions are under changes in this input power 
spectrum. We have repeated our analysis using input power spectra with primordial power 
spectra having slopes n = 1.5 and n = 0.5 instead of the standard value n = 1. To give 
an idea of the magnitude of this variation in n, we remind the reader that the COBE data 
determine n with a one-sigma uncertainty of about 0.3 (Tegmark & Bunn 1995). The amplitude 
of the predicted structure increases with increasing n, as one would expect. When we vary n by 
±0.5, the two highest peaks in the lowest panel of Figure 7 change by ±11 /iK and ±9 /iK, or 
approximately 15-20%. The location of the peaks does not change significantly. The presence 
and approximate height of these two peaks can therefore be regarded as robust predictions. 

5 Conclusions 

Wiener filtering is a promising tool for the analysis of CMB sky maps. Filtering provides a 
significant improvement in the signal-to-noise ratio in the regions covered by the raw data, and 
allows some small amount of information to be reconstructed about the anisotropy within the 
Galactic cut region. As Figure 2 shows, the most prominent hot and cold spots in the filtered 
map carry very high statistical significance. 

This reduction in noise is not without a price. In order to apply the Wiener filter, one needs 
to assume a power spectrum, and the cleaned data therefore depend on more assumptions than 
do the raw data. However, the large-scale features in the filtered map are only weakly dependent 
on the choice of power spectrum. 

One great advantage of the Wiener-filtered map is that the statistical properties of the 
residuals are known precisely, assuming the underlying model is correct. It is therefore possible 
to assess the goodness of fit of any other data set with the combined DMR data and assumed 
theory. In addition, it is easy to make constrained realizations of the microwave sky with the 
correct a posteriori probability distribution. This technique can in principle be used to perform 
Monte-Carlo simulations of the CMB for comparison with experiments on similar angular scales. 

The Tenerife experiment seems to be the most promising experiment on which to test the 
possibility of using the Wiener filtered maps to make predictions. Unfortunately, the triple-beam 
strategy used in the experiment drastically reduces Tenerife's sensitivity to the modes that are 
constrained well by the DMR. It is possible, however, to predict some fraction of the Tenerife 
signal from the filtered DMR data. By incorporating this information into a Tenerife analysis, 
it should be possible to place somewhat tighter constraints on models than would be possible 
with the Tenerife data alone. 

In particular, the DMR data predict two strong features in the Tenerife data at right ascen- 
sions around 220° and 255° and declinations around 35° (see Figures 7 and 8). The statistical 
significance of these features is high enough that if they are not found to be present in the real 
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data, we will have strong reason to conclude that something is wrong with either the assump- 
tions underlying the prediction or the data from one of the experiments. Specifically, if the 
Tenerife experiment were to fail to see these features, we would be led to one of the following 
conclusions: 

1. The CMB anisotropy is not well described by Gaussian statistics. (The Wiener filtered 
map does not depend on the assumption of Gaussian statistics, but the uncertainties in it 
do.) 

2. The CMB anisotropy has a power spectrum drastically different from that usually assumed. 

3. One of the two experiments is contaminated by systematic error or foregrounds. 

The predicted features are fairly robust to reasonable changes in the CMB power spectrum 
on COBE scales, so the only way the second option could be correct is if the power spectrum 
has a drastic rise at high enough / to be invisible to COBE. It is somewhat difficult to see how 
the third option could be correct: the COBE data have been extensively checked for systematic 
errors and foregrounds (Bennett et al. 1992, Kogut et al. 1992), and it would be difficult to 
explain the failure of Tenerife to see the predicted signal by assuming a contamination of the 
Tenerife data: contaminants are much more likely to induce spurious signals than to mask real 
ones. Naturally, any of these three conclusions would be of enormous interest. 

We would like to thank M. White for supplying us with the CDM angular power spectrum. 
Y.H. acknowledges the hospitality of the Center for Particle Astrophysics. Y.H. has been sup- 
ported in part by The Hebrew University Internal Funds (grant 53/94) and by the Israel Science 
Foundation administered by the Israeli Academy of Sciences and Humanities (grant 590/94). 
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Figure Captions 

Figure 1: The Wiener-filtered DMR data 

The result of applying a Wiener filter to the two-year COBE DMR data. The input 
power spectrum is that of a standard cold dark matter model. The contour levels are 
at ±25, ±50, ±75 /iK, with negative values indicated by dashed lines. The region of 
sky observed by the Tenerife experiment, i.e., the region with equatorial coordinates 
140° < R.A. < 260° and 35° < 8 < 45°, is also marked on the greyscale map. In 
Section 4 we will use the filtered DMR map to make predictions for this experiment. 
The part of the strip near 260° in right ascension is split in half in this projection. 



Figure 2: Signal-to-noise ratio in the filtered map 

The ratio of the filtered map in Figure 1 to the uncertainty is plotted. The contours 
take values ±2, ±3, ±4, with dashed lines representing negative values. 



Figure 3: Constrained realizations 

Four constrained realizations of the microwave sky are shown. The realizations 
represent typical expected fluctuations about the Wiener-filtered map in Figure 1, 
and were drawn from the Gaussian probability distribution described in the previous 
section. Each realization is shown in both greyscale and contour maps. The contour 
levels are ±50, ±100, ±150 /iK, with negative values indicated by dashed lines. 
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Figure 4: Power spectrum of the Wiener-filtered map 

The power spectrum A\ defined in equation (22) is shown. The standard CDM 
input power spectrum is also plotted for comparison. 



Figure 5: DMR and Tenerife window functions. 

The power-spectrum window functions Wi are shown for both the COBE and Tener- 
ife experiments. The one with little power at low / is Tenerife. The dashed line shows 
/(/ + 1)C/, where C\ is an arbitrarily-normalized standard cold dark matter power 
spectrum. 



Figure 6: Tenerife predictions 

A map of the most probable (a posteriori) Tenerife signal over the region 120° < 
RA < 260°, 35° < S < 45°. The greyscale ranges from -60 /iK to 60 fiK, and the 
uncertainty associated with each point is approximately 28 jiK. 



Figure 7: Tenerife predictions with uncertainties 

The solid lines show slices through the map in Figure 6, and the dashed lines show 
plus and minus one-sigma errors. The top, middle, and bottom panels correspond to 
declinations of 45°, 40°, and 35° respectively. 



Figure 8: Constrained realizations of the Tenerife data 

We show five constrained realizations of the Tenerife predictions to assist the reader 
in qualitatively assessing the uncertainties in the predictions in Figures 6 and 7. The 
top panel shows five realizations covering the same range as Figure 6. The greyscale 
in this map ranges from -110 to 100 /iK. The bottom panel is a strip through the 
five realizations at declination 35°, and is thus to be compared with the lowest panel 
of Figure 7. 
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